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Abstract 

In this paper we investigate the short-time decoher- 
ence from Ohmic and 1/f noise of single Josephson 
charge qubit (JCQ). At first, we use the short-time 
approximation to obtain the dynamics of the open 
JCQ. Then we calculate the decoherence the mea- 
sure of which is chosen as the maximum norm of the 
deviation density operator. It is shown that the de- 
coherence from 1/f noise plays the central role. The 
total decoherence from Ohmic and 1/f noise is se- 
rious at present experiential conditions according to 
the DiVincenzo criterion. 
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1 Introduction 

Quantum bit (qubit) is a key block for building quan- 
tum computers. Among various realizations of the 
qubit for quantum computations that based on the 
Josephson junctions are considered to be particu- 
larly promising candidates because of their scalabil- 
ity, established fabrication techniques, and flexibil- 
ity in designs. Many kinds of the superconducting 
Josephson-junction-qubit models are proposed in last 
years. They are Josephson charge [Q, flux PJ, phase 
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and hybridized 0] qubits. According to DiVin- 
cenzo [5] proposal, the qubit for quantum computa- 
tion must satisfy five criteria one of which is the low 
decoherence criterion. An approximate benchmark 
of the criterion is a fidelity loss no more than ~ 10~ 4 
per elementary quantum gate operation. The deco- 
herence of the Josephson single qubit and coupling 
qubits has been investigated widely in last years |ü]. 
In these researches the dynamics of the qubits inter- 
acting with their environment has been treated by the 
perturbation method, the path integral method [Zj, 
and others [Hj- In the researches in general some ap- 
proximation scheme must be used. The most familiar 
and frequently used one is the Markov approximation 

m 

It has been pointed that the Markov approxima- 
tion can not be used in low temperature and for short 
cycle times of quantum computation [TJJ. However, 
the qubits based on the Josephson superconducting 
junction should be worked in low temperature, such 
as few 10 mK. So the Markov approximation can 
not be used in the investigations of decoherence for 
Josephson qubits. Fortunately, a short-time approxi- 
mation scheme m] being fit for the investigations has 
been developed recently by Privman et al. ^| Ej ■ 
The approaches of V. Privman and his co-worker are 
rather formal and universal. It is also very interest- 
ing to investigate the decoherence of a concrete and 
physical qubit model because it can help us to know 
that whether the model satisfy the DiVincenzo low 
decoherence criterion or not. If the model can not 
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satisfy the criterion, the qubit has then not qualifica- 
tion for being a quantum computation hardware. 

In this paper we shall investigate that if the JCQ 
model satisfy the DiVincenzo low decoherence cri- 
terion. The same problem has been investigated in 
our previous works ^] where we suppose that 
there is only the Ohmic noise in the environment of 
JCQ. However, recent several experiments ( see [Tïï] 
and within ) with Joscphson junction circuits have 
revealed at low freqüències the presence of 1// noise 
and it is shown that the 1// noise plays a major role 
in destructing the coherence of qubit in solid-state 
systems. In this paper, we shall carefully estimate the 
decoherence of the JCQ in Ohmic and 1// noise envi- 
ronment. Here, the measure of decoherence is based 
on a Standard operator norm \\A\\ in the theory of 
linear operators |17| . The short-time approximation 
of the split-operator will be used in our investigation. 

2 Open JCQ model 

The single JCQ Hamiltonian is [T%1ITÜ] 



is ESj then 



H s = AE C (n — n g ) — Ej cos (p. 



(1) 



Here, E c — e 2 /2 (C g + Cj) is the charging energy; 
Ej = I c h/2e is the Josephson coupling energy |2U) . 
where I c is the critical current of the Josephson junc- 
tion, h the Planck's constant divided 2ir, and e the 
charge of electron; n g = C g V g /2e = Q g /2e is the di- 
mensionless gate charge, where C g is the gate capac- 
itance, V g the controllable gate voltage. The numbcr 
operator n of (excess) Cooper pair on the island, and 
the phase operator ip of the superconducting order 
parameter, are quantum mechanically conjugate 
The unavoidable noise of the environment may lead 
to the dissipation. In general, one takes into account 
that the environment itself is also a quantum system 
with a large number degrees of freedom. Usually it 
is modeled by a bath with a large set of harmònic 
oscillators J2HI35] 



(2) 



each of which interacts weakly with the system of 
interest. The whole Hamiltonian of the system-bath 



H 



2m n 



1 

+ -m a uj^ I x a 



j-2en— 



(3) 



Here, the coupling operator of the environment is 
X = ^2X a x a . If the Josephson coupling energy Ej 
is much smaller than the charging energy E c , and 
both of them are much smaller than the supercon- 
ducting energy gap A, the Hamiltonian H s of Joseph- 
son junction can be parameterized by the numbcr 
of the Cooper pairs n on the island. When the 
temperature T is low enough the system can be re- 
duced to a two-state system (qubit) because all other 
charge states have much higher energy and can be 
neglected. So the Hamiltonian of the system can ap- 
proximately reads H s = —\B z a z — ^B x a x , where 
B z = E c h (1 — 2n g ) and B x = Ej. If we choose the 
working point to make n g half-integer, say n g = 1/2, 
we can obtain that 



Here, 



H^Hs + Hj + Hb- 



= --B x a x , H B = J2 a tw a ala a , 
= -a z VhY^,9a («4, + a a ) , 



(4) 



where 



m a uj a \ Cj 



(5) 



(G) 



3 Ohmic and 1/f noise 

Recent studies, both the experimental and theoreti- 
cal, suggest that the serious noise sources in Joseph- 
son devices include two kinds of fluctuations one of 
which with ohmic spectrum the other with 1/f spec- 
trum |16j . In the following, we introduce the spectral 
dcnsities of the Ohmic and 1/f noise. According to 
|21j , from Eq. J3J , we can express the spectral density 
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of the noise as the JCQ circuit the typical impedance of the control 

J g 



2 2 line is R ps 50 Q, and C g ps 1CT 18 F, Cj ps 10~ 16 F, 



J x (u) = E ( 2e-§-") i^L^ ( w - Wq ) so we can obtain 77 ps 10 6 

2 q V w / rn a uj a The 1// noise is considered deriving from the back- 

_ —hJ2 \g a \ 2 S (uj — Lü a ) ground charge fluctuations in the circuits. It can be 

2 a expressed an effective noise of the gate charge, i. e., 

7r^ í°° )\ 2 5 ( ) ^ = Oífe /lj. Recent experiments proposed at 

2 1 c /_ 00 a W relevant temperatures a/ ~ 10~ 7 — 10~ 6 . According 

7r 2 to Ref. 23 this noise can be translated into the fluc- 

= -KD(u>)\g{u>)\ , (7) tuations of X, namely, 

where D (u) is the density of the states of environ- g x = Ej/uj, (13) 
ment modes. On the other hand, the spectral density 

is related to the power spectrum as [21ESI where Ef = 4E c ^aJ. Comparing the Eq. I|13(l with 

Eq. JSJ we have 



S x H = J x M S coth (w/3/2) . (8) 



l&E^a 



...... . "-""/ 

Here, /3 = h/k^T, where kn is the Boltzmann con- Jx (w)| f = -7-7 — . , . (14) 

X11 11 1 1 • 1. m ftw coth (cj/í/2) 
stant. It has bccn strcsscd that the relationship J8J) 

are satisfied not only for the linear coupling model So 

but also for the nonlinear onc. To the Ohmic noise 

case, from the fiuctuation-dissipation theorem one n r. ,\ „2 / \ 1 32E c af 

has ESI K >9 { >lf nh*u coth (w/3/2) 



_ / 

ç r M f 2eC *V ç , M a; coth M/2)' 
„ where 

2 «/I p2 



(15) 



= f^rO HeZt(w)w»coth(w/8/2). k = ps 1.5 x 10 25 . (16) 

\ Cj J hn 

(9) In fact, in a JCQ circuit these two kinds of noise 

are existed at the same time. So the spectral density 
Then, setting Re Z t (w) ps ü, one has rf thg total noise . g 



J x (uj)\ = (2e) 2 (^) Ru. (10) J(w)= Jx(w)| + JxH| r (17) 

__ So we have 
So by using Eq. (Jjj we have 

,„ ,2 /r X 2 D(u)g 2 (u) = D(u)g 2 (u>)\ + D(u>)g 2 (u)\ f . 

D(u)g 2 (u;)\ = 4 { -f- flw. (11) (18) 

V J / In the following, we shall use above three kinds of 

Considering the cutoff of the frequency w we can set different noise of the environmcnt invcstigating the 

decoherence of the JCQ. 



D(u)g i (u)\ =riuexpl — 1 . (12) 

4 Short-time dynamics of JCQ 

Here, lü c is the cutoff frequency and ri — 4-^- ( Sr- ) , D r - 1 • ,1 1 1 un- 1 • 1 

' H J ' «q \w/ Reforç studymg the decoherence we shall mvestigate 

where Rq — h/ (2e) 2 ps 6.5 kfí |24| . In the model of the dynamics of the open JCQ with the short-time 
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approximation. Suppose the initial state of the JCQ- 
bath be -R(O) = p(0) £5 G, where p(0) is the initial 
state of JCQ and is the initial state of the envi- 
ronment. We set O is the product of the bath modes 
density matrices 9k- In the initial states, each bath 
mode k is assumed to be thermalized, namely, 



-/3M k 



'fc = 



(19) 



where Mj. = u)ka\a,k- 
JCQ-bath is then 



Tr fe (e-P M >y 
The evolution operator of the 



U 



^-íHt/K 



^-iiH.+Hr+Hn^/H 



(20) 



Due to non-conservation of H s in this system, the 
evolution operator cannot be in a general way ex- 
pressed as e -iH B r/n. e -i{H I +H B )T/h_ But m the sort _ 

time approximation, the operator can be approxi- 
mately expressed as j^S] 



U 



e -iH B r/2h e - 



■i(Bi+H B )r/h e -iH s r/2h + , T 3y ^1) 



It has been proved that the expression is accurate 
enough for the time being short to the characteristic 
time. In the following we only investigate the case 
that the system evolute within the time r < 10 ps, 
the characteristic time of the single JCQ is r = 12.7 
ps (when Ej — 51.8 peV) |15| . So the elements of the 
reduced density matrix p (r) in the basis of operator 
H s can be expressed as 



where p m = p 00 (0) , p n = p n (0), p w = p w (0) , and 



B 2 (í) 



C(t) 



fe 



2 U k t @UJ k 

i coth . 



E 

fe 



\9k\ 



(ujkt - sinwfcí) . 



(24) 



When the summation in Eq. (|24|l is converted to 
integration in the limit of infinite number of the bath 
modes, one has 



N 2 u 2 sm' 



2WÍ coth^ 



B 2 (t) = 8 J ^M 3 H 2 .. .... ... 

(25) 

for the real g (u>) . In fact this B 2 (t) is the mean- 
squared value of the magnitude of the phase noise 
in Ref. From the deriving of B 2 (t) we know 

that it results from the choice of bath model and the 
coupling forms of JCQ-bath rather than the choice 
of the approximation scheme. But the evolution of 
the matrix elements as Eqs. I)23|l is resulted from the 
short-time approximation and it can not be obtained 
by using the Markov approximation and others. By 
using above results we can yield a good quantitative 
estimation to the decoherence behavior of JCQ in a 
short time. 



5 Decoherence of JCQ 



It is shown that the measure 



Pn 



TÏB {<p v 



■/2h e -i{H,+H E 



r / n e- iH > T/2h R(0) 



£>(í) = su P (|Kt,p(o))y 

P(0) 



(26) 



e iH s T/2h e HH I +H B )T/h e iH B r/2h 



(22) 



where {m, n} = or 1. The p is a matrix with 2 by 2. 
In the following we set t — t/H. Through some cal- 
culations we can obtain the evolution of the density 
matrix elements p 10 (t) and p xl (t) as |14l 115) 



is suitable for estimating the decoherence of the qubit 
gates Here, the norm ||tr|| A is defmed as 

i 

'(tp\<r\tp) 



= sup 



(27) 



For a qubit, it is 



P W (*) = 2P10 ( 1 - e - S2(t) + ^ + 



ctio | 2 + klll • 



( o 11 (í) = ip Q0 ^1 — e~ s2 ^^j + TjPn + e fi2< ' í ^ , Here, the deviation operator a is defined as 

(23) (T(r)=p(r) _ p i (r)) 



(28) 
(29) 
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where p (r) and p l (r) are density matrixes of the 
"real" evolution (with interaction) and the "ideal" 
one (without interaction) of the investigated system. 
The evolution of the closed JCQ is p\ x (t) — p 11: and 
Pio (*) = Pio e1,tEj ■ So for the open JCQ we have 



(30) 



«no (t) = 5P10 (l - e- fl 'W) (1 ^ e*^ 
^11 (*) = ~ (l - e _S2(í) ) (Poo ~ • 



Thus, we have 



kWII 



l('--* w ) 



2 

Ooo 



Pn) 2 + 4|p 10 | 2 sin 2 ^} T 



From Eq. (|31[) we know that a pure state 



Pi (0) 



1 




or p 2 (0) = 




1 



(31) 



(32) 



will make the ||cr(í)|k into D (i). In the following 
calculations we use the initial state p 1 (0) and choose 
the Josephson energy Ej — 51.8 /icV according to 
Ref. p. 



Fig.l 



Fig.2 



Fig.3 



noise. Here, we confine the freqüències of bath modes 
in 1 GHz^ 50 GHz this is a very broad band of 
freqüències of the environment noise. The lowest 
line in Fig.l plots the decoherence at temperature 
T x = 0.0300 K. It is shown that when t < 0.1, 
namely, r < 6.582 x 10~ n s the decoherence is much 
smaller than DiVincenzo criterion 10 -4 . The middle 
and the upmost lines plot the decoherence at tem- 
perature T 2 = 0.1500 K and T 3 = 0.1875 K. From 
the plots we see that the decoherence will increase 
with the increasing of the temperature T but it can- 
not increase to the maximum valuc of the DiVin- 
cenzo criterion permitting at T = 0.1875 K within 
t = 0.02 ( r ps 12.7 ps). Fig.2 shows the decoher- 
ence from 1/f noise of the JCQ in different a/. It is 
shown that this decoherence mainly derived from the 
low frequeney parts of the 1/f noise and it is sensi- 
tive to the changing of otf. In the calculations, we 
take the freqüències of the environment modes from 
1 kHz to 1 GHz the numerical simulation shows that 
the spectrum is broad enough. In the Fig.2 the low- 
est, middle, and upmost lines plot the decoherence at 
otf = 1.0 X 10~ 7 , 1.1 x 10" 7 and 1.2 x 10~ 7 . A further 
calculation shows that if only a/ > 5.0x 10 ~ 8 the de- 
coherence will affect the JCQ as the qubit for making 
quantum computers because the decoherence is big- 
ger than the DiVincenzo low decoherence criterion 
in the single operation time. So we should decrease 
the a/ in order to make JCQ be a qubit for quan- 
tum computation. In Fig.3 we simulate three kinds 
of cases where both the Ohmic and 1/f noise are con- 
sidered. Here, we set T = 30 mK and a/ = 5 x 10~ 8 
(the lowest line), a/ = 4 x 10~ 8 , (the middle line) 
and otj = 3 x 10~ 8 (the upmost line). The freqüèn- 
cies of bath modes is set in 1 kHz^ 1 GHz for 1/f 
noise and in 1 GHz^ 50 GHz for Ohmic noise. It 
is shown that when the otj do not exceed 5 x 10~ 8 , 
within t = 12.7 ps the decoherence from the Ohmic 
noise and 1/f noise is endurable. 



By using Eqs. I|11I15I18|I . and (j2U we can 

plot the decoherence of JCQ versus time t. In Fig.l 
we plot the decoherence from Ohmic noise of the 
JCQ versus time t. This decoherence is mainly de- 
rived from the higher frequeney parts of the Ohmic 



6 Conclusions 

In this paper we investigated the short-time deco- 
herence results from the Ohmic and 1/f noise of the 
JCQ. It is shown that the decoherence from 1/f noise 
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is larger than that from Ohmic noise to the JCQ 
model. To the Ohmic noise the higher frequency 
parts play a major role to the decoherence. We take 
a larger range of the freqüències of the environment 
modes in our numerical simulation than usually pro- 
posed. It is shown that in usually experimental tem- 
peratura this decoherence is not serious comparing 
to the DiVincenzo criterion. Unlikc the Ohmic noise 
case, to the 1// noise the decoherence is mainly de- 
tcrmined by the lower frequency parts. It is shown 
that when aj < 5 x 10~ 8 , within r = 12.7 ps this de- 
coherence is also not serious according to the DiVin- 
cenzo criterion. However, the present experiments 
show that the value of a/ is between 10~ 7 and 10 -6 . 
Thus, if we wish to make the JCQ become a optimal 
qubit model we should decrease the value of a/ . The 
value of a/ is related to the temperature. It may 
decrease with the decreasing of temperature. Care- 
fully finding out the correlation of the value of ctf 
and other parameters (include the temperature) is a 
important task for our investigating the decoherence 
of JCQ model in the future. 
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7 Figs. captions 

Fig.l The decoherence from Ohmic noise within short 
time. The lowest, middle and upmost lines cor- 
respond to T 2 = 0.0300 K, T 2 = 0.1500 K and 
T3 = 0.1875 K. Here, the freqüències of bath modes 
is set in 1 GHz- 50 GHz. 

Fig.2 The decoherence from 1// noise within short 
time. The lowest, middle and upmost lines corre- 
spond to otf — 1.0 x 10~ 7 , a/ = 1.1 x 10~ 7 , and 
ctf = 1.3 x 10~ 7 . Here, the freqüències of bath modes 
is set in 1 kHz~ 1 GHz. 

Fig.3 The decoherence from Ohmic noise and 1// 
noise within short time. The lowest, middle and 
upmost lines correspond to a/ = 3 x 10~ 8 , o.f — 
4 x 10~ 8 , and a/ = 5x 10~ 8 . Here, the temperature 
is set T = 30 mK, the freqüències of bath modes is 



set in 1 kHz^ 1 GHz for 1// noise and in 1 GHz~ 50 
GHz for Ohmic noise. 
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